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RESONANCE PARAMAGNETIC RELAXATION AND ALIGNMENT OF SMALL GRAINS 
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ABSTRACT 

We show that the energy-level splitting arising from grain rotation ensures that paramagnetic dissipa- 
tion acts at its maximum rate, i.e., the conditions for paramagnetic resonance are automatically fulfilled. 
We refer to this process as "resonance relaxation" . The differences between the predictions of classical 
Davis-Greenstein relaxation and resonance relaxation are most pronounced for grains rotating faster 
than 1 GHz, i.e., in the domain where classical paramagnetic relaxation is suppressed. This mechanism 
can partially align even very small grains, resulting in linearly polarized microwave emission which could 
interfere with efforts to measure the polarization of the cosmic microwave background. 

Subject headings: ISM: Atomic Processes, Dust, Polarization; Cosmic Microwave Background 



1. INTRODUCTION 

Experiments to study the cosmic background radiation 
have stimulated renewed interest in diffuse galactic emis- 
sion. Recent maps of the microwave sky brightness have 
revealed a component of the 10-100 GHz microwave con- 
tinuum which is correlated with 100 thermal emission 
from interstellar dust (see review by Draine & Lazarian 
1999a). Draine & Lazarian (1998a, b, henceforth DL98a,b) 
attributed this emission to electric dipole radiation from 
small (< 10~ 7 cm) rapidly rotating grains. Recent ob- 
servations by de Oliveira-Costa et al. (1999) support this 
interpretation. The question now is whether these small 
grains are aligned and their emission polarized. 

One process that might produce alignment of the ul- 
trasmall grains is the paramagnetic dissipation mecha- 
nism suggested by Davis and Greenstein (1951) to explain 
the polarization of starlight. The Davis-Greenstein mech- 
anism is straightforward: the component of interstellar 
magnetic field perpendicular to the grain angular veloc- 
ity varies in grain coordinates, resulting in time-dependent 
magnetization, energy dissipation, and a torque acting on 
the grain. As a result grains tend to rotate with angular 
momenta parallel to the interstellar magnetic field. 

Although recent research (Draine & Weingartner 1996, 
1997, Lazarian & Draine 1999a, b) suggests that paramag- 
netic alignment may not be the dominant alignment mech- 
anism for a ^ 10~ 5 cm grains, it may be effective for small 
(a ^ 5 x 10 -6 cm) grains. 

In the present paper we claim that the traditional pic- 
ture of paramagnetic relaxation is incomplete, as it dis- 
regards the splitting of energy levels that arises within a 
rotating body. Unpaired electrons having spin parallel and 
antiparallel to the grain angular velocity have different en- 
ergies resulting in the Barnett effect (Landau & Lifshitz 
1960) - the spontaneous magnetization of a paramagnetic 
body rotating in field-free space. Therefore the implicit as- 
sumption in Davis & Greenstein (1951) - that the magne- 
tization within a rotating grain in a static magnetic field is 
equivalent to the magnetization within a stationary grain 
in a rotating magnetic field - is clearly not exact. 

In what follows we show that a very important effect 



due to rotation has thus far been overlooked. This ef- 
fect, which we term "resonance relaxation" , leads to en- 
ergy dissipation - and grain alignment - which is much 
more rapid than the classical Davis-Greenstein estimate 
when the grain rotates very rapidly. 

2. DAVIS-GREENSTEIN THEORY 

Paramagnetic dissipation in a stationary grain depends 
upon the imaginary part of the magnetic susceptibility 
X", which characterizes the phase delay between the grain 
magnetization and the rotating magnetic field. Due to 
this delay a grain rotating in a static magnetic field expe- 
riences a decelerating torque: the energy dissipated in the 
grain comes from rotational kinetic energy. The Davis- 
Greenstein alignment time scale is 



tug ~ 3 x 10 yr 
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where K(uS) = x"(lu)/u) (see Davis & Greenstein 1951) 
Following DL99b we estimate 
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Very small grains are expected to be paramagnetic both 
due to presence of free radicals, paramagnetic carbon rings 
(see Altshulcr & Kozyrev 1964) and captured ions. The 
spin-spin coupling time 

h 1Q _ 9 / l0 21 cm- 3 \ g 

3.8n P 3/i| \ n p J 

(see DL98b) where /is is the Bohr magneton, and n v ~ 
10 21 cm -3 is the concentration of unpaired electrons, 
greater than in coals (Tsvetkov et al 1993), but less than 
the concentration of free radicals envisaged by Green- 
berg (1982). Eq. (3) then predicts a cut-off frequency 
t'cut = (27TT2) -1 w 0.1 GHz. In the extreme case where 
~10% of the atoms are paramagnetic one can get is cut as 
large as 1 GHz but hardly any higher. 

3. BARNETT EFFECT AND BLOCH EQUATIONS 

The Barnett effect states that a body rotating with ve- 
locity u> develops a magnetization 

M = - X o— =XoH BE ■ (4) 
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where xo is the static susceptibility, and Hbe is the 
"Barnett-equivalent" field. The essence of the Barnett ef- 
fect is easily understood: a rotating body can decrease its 
energy, while keeping its angular momentum constant, if 
some of the angular momentum is taken up by its unpaired 
spins. By flipping one spin of angular momentum H/2, the 
system can reduce its rotational kinetic energy by Hu. 

Although the Barnett effect has been long known in 
physics (see Landau & Lifshitz 1960), its importance in the 
context of interstellar grains was only appreciated recently 
(Dolginov & Mytrophanov 1975; Purcell 1979; Lazarian 
& Roberge 1997; Lazarian & Draine 1997, 1999a). In the 
present paper we discuss a hitherto- unrecognized aspect of 
the Barnett effect, namely its influence on paramagnetic 
dissipation in a rapidly rotating grain. 

Rotation removes the spin degeneracy of the electron 
energy levels. The energy difference between electron spin 
parallel or antiparallcl to tu provides a level splitting cor- 
responding to hu = Insofar as the energy levels 
and magnetization are concerned, rotation of the grain is 
analogous to application of the "Barnett equivalent" field. 

Now consider a (weak) static magnetic field Hi at an 
angle 9 to u>. In grain coordinates, this appears like a 
static field Hi cos 9 plus a field Hi sin 9 rotating with fre- 
quency u. This rotating field can be resonantly absorbed, 
since the energy level splitting is exactly hu. 

The Bloch equations (Bloch 1946) are useful for describ- 
ing both resonant and nonrcsonant absorption (see Pake 
1962). These phcnomenological equations reflect the ten- 
dency of the magnetization M to precess and to tend ex- 
ponentially towards its thermal equilibrium value Mq. In 
a stationary grain Mo = XoH, where \o is the static para- 
magnetic susceptibility and H is the external field. In the 
case of a rotating grain the magnetization arises due to 
the Barnett effect and is therefore given by eq. (4). 

In what follows we assume that the grain magnetiza- 
tion is directed along the z-axis and is dominated by the 
Barnett effect (the magnetic field for a grain rotating at 
20 GHz is 7.2kG, much greater than the internal field in a 
paramagnetic grain). Changes in M along this axis involve 
changes in the energy of the spin system, thus require spin- 
lattice interactions, and therefore occur on the spin-lattice 
relaxation time scale t\ (Athcrton 1973). Changes in M 
in the perpendicular direction only slightly perturb its di- 
rection but not its magnitude. The interactions within 
the electron spin system are sufficient to deflect the direc- 
tion of magnetization and these perturbations relax on the 
spin-spin relaxation time t 2 . 

Using _L for the x and y components, the Bloch equa- 
tions in the presence of the interstellar magnetic field Hi 
are (see Morrish 1980) 

M ^ (5) 
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where dM±/dt represents the motion of in body co- 
ordinates and 7 = e/2m e c = 8.8 x 10 6 s _1 G _1 . 

Consider a reference frame (i,j,k) rotating with the 
grain at angular velocity u. In this frame the stationary 
interstellar magnetic field is 

Hi = [iHi cos(wt) - ]H x sin(wt)] sin 6 + kHi cos 9 . (7) 



As the magnetization in the z-direction xoHi cos 9 is much 
smaller than xo^be, it is disregarded in our treatment. 

For M z the stationary solution is M z = — xoHbe/ fst 
where / st = l+"f 2 g 2 TiT 2 H 2 sin 2 9 1 while M x and M y oscil- 
late with a 7r/2 lag with respect to the interstellar magnetic 
field. For instance, M x = xout 2 Hi sin(u;t)// st . Therefore 



X 



Xo 



UT 2 



l+l 2 g 2 s TiT 2 Hlsm 2 9 



(8) 



which coincides with the expression for \" for electron 
paramagnetic resonance in a stationary sample when the 
frequency of the oscillating field Hi is equal to the res- 
onance frequency. In our problem the only relevant fre- 
quency is the frequency of grain rotation. Therefore it is 
not accidental that the paramagnetic relaxation is "res- 
onant" when the grain rotates in the external magnetic 
field. The term 7 s griT 2 H 2 sin 2 9 in the denominator al- 
lows for "saturation" when the energy dissipated in the 
spin system raises its temperature due to slow spin-lattice 
coupling; below we show that this term can be important 
for very small grains, for which n can be large. 

An important difference between paramagnetic reso- 
nance in an external magnetic field and the resonance re- 
laxation discussed here is that the "Barnett equivalent" 
magnetic field is different for species with different mag- 
netic moments. Therefore, unlike paramagnetic resonance, 
resonance relaxation happens simultaneously to species 
with completely different magnetic moments and g-factors. 
For example, the conditions for electron spin resonance 
and nuclear magnetic resonance are satisfied simultane- 
ously when a grain rotates in a static weak magnetic field. 

4. SPIN-LATTICE RELAXATION 

For a spin to flip in a rotating grain, it is necessary for 
the total energy in lattice vibrations to change by hu. Be- 
cause the density of states is finite, it may not be possible 
for the lattice vibrations to absorb an energy hu. 

To estimate the lowest vibrational frequency w m i n we 
note that the lowest frequency bending mode of the 
coronene molecule (C24H12) is u m i n = 1.9 x 10 13 s _1 
(Cyvin 1982). Coronene has an effective radius a — 4.3 A 
for an assumed density p — 1.5 gem -3 ; if u m - m oc a -1 , then 
w m in = 8.3 x 10 12 oC 7 s _1 , large compared to kT/h. Thus 
for ultrasmall grains one cannot use spin-lattice relaxation 
times n measured for macroscopic samples. 

We obtain an upper estimate for ti by appealing to 
the Raman scattering of phonons (see Waller 1932, Pake 
1962): annihilation of a vibrational quantum hu' and cre- 
ation of a quantum h(u — u'). Elastic lattice vibrations 
of all frequencies participate and therefore the relaxation 
is present for small grains. After integrating over vari- 
ous modes of vibrations one gets the probability of the 
spin-lattice transition via Raman scattering (Al'tshuler & 
Kozyrev 1964) 



rf 1 « K 2 (kT/h) m+1 J m /pcf , (9) 

where p is the density, c s the sound speed, K 2 is a function 
that depends on the density of states, and J m (m =6 or 
8) is an integral over the body's phonon frequencies: 
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where 9d is the Debye temperature. The conventional 
treatment assumes the body to be infinite with the inte- 
gration extending down to Ti = 0. In this case, and for 
T <c ^o/m, we have J m w m!£(m), where £ is the Rie- 
mann zeta function. For a grain of size a, 

7ia; m i n 63K 
J; = — t— « • (10 

k a_7 

For T <C Ti/m we have </ m ~ (Ti/T) m exp(— Tj/T). The 
ratio of the Raman spin-lattice relaxation in a small grain 
at temperature T <^Ti jm to such relaxation in an infinite 
body at 77 'K £ On/m is then 

„ (™)™ +I (|)™ xp(r , /T)m , ((ra) . 

Data in Al'tshuler & Kozyrev (1964) suggests that ionic 
crystals have a spin- lattice relaxation time Ti !0O (77K) « 
10 -6 s. If due to the Raman process, then we would es- 
timate that a macroscopic sample at T = 4K would have 
n rs 28 s if m = 6, or 2 x 10 4 s if m = 8. ^From eq. 
(11), a grain with a = 10~ 7 cm at T = 4K would have 
n w 3.1 x 10 5 s for m = 6, or 2.7 x 10 7 s for m = 8. Grains 
with a < 10~ 7 cm would have even larger values of t\. 
If we adopt t-i from eq. (3), then 

' Hi \ 2 sin 2 fl 

10 6 s J \5^GJ 2/3 

so we see from (8) that saturation may be important for 
a < 1CT 7 grains even in the ~ 5[iG fields in diffuse inter- 
stellar gas. 

How reliable is our above estimate for n? Our calcu- 
lations were based on so-called Waller theory, which fre- 
quently overestimates T\ by a factor up to 10 8 (Pake 1962). 
If the dependence of the spin-lattice relaxation time on 
temperature is different from that given by Eq. (9) our es- 
timates of T\ at 4K would be very different. Ultimately we 
require laboratory measurements of n in small particles of 
appropriate composition. 

5. GRAIN ALIGNMENT 

Paramagnetic alignment of grains with a given axis ra- 
tio depends on two parameters: the ratio Td/T rot of grain 
vibrational and rotational temperatures, and the ratio of 
the alignment time tdg to the rotational damping time fa- 
The rotational damping time t^ depends on various pro- 
cesses of damping and excitation (e.g. collisions with ions 
and neutrals, plasma drag, and emission of photons) dis- 
cussed in DL98b. Assuming that the paramagnetic torque 
only marginally reduces u, we follow the analysis of DL98b 
to obtain w as a function of a. 

In our calculations we assume that a grain spends most 
of its time between thermal spikes with a vibrational tem- 
perature T Rj 4K (cf. Rouan et al. 1992), since the time 
between photon absorptions is ~ 10 9 s, while the grain 
cools much more rapidly. Note that photons usually con- 
tribute marginally to the disorientation of the grain angu- 
lar momentum J (see Fig. 5 in DL98b). 

Let 9 be the angle between J and the interstellar mag- 
netic field B. Fig. 1 presents the measure of alignment 
a = | (cos 2 9—1/3), for grains in the cold neutral medium. 



For our estimate we used standard formulae for paramag- 
netic alignment of angular momentum with B (Lazarian 
1997, Roberge & Lazarian 1999), which for weak alignment 
provide a ~ 2/15[l — (1 +ri)/(l +r)], where r = tdgAcI 
and t = T/T rot . 

The discontinuity at « 6 x 1CP 8 cm is due to the assump- 
tion that smaller grains are planar, and larger grains are 
spherical. The degree of polarization of rotational electric 
dipolc emission p « a cos 2 ip, where ip is the angle be- 
tween B and the plane of the sky. pj cos 2 ip as a function 
of frequency is also shown in Fig. 1. The dipole rota- 
tional emission predicted in DL98a,b is sufficiently strong 
that polarization of a few percent may interfere with ef- 
forts to measure the polarization of the cosmic microwave 
background radiation. It worth noting that the degree of 
microwave polarization is sensitive to the magnetic field 
intensity (through tug)- 

6. DISCUSSION 

We have discussed a new gyromagnetic effect - reso- 
nance relaxation - which is closely related to normal para- 
magnetic resonance, and arises naturally whenever a body 
rotates in a weak magnetic field. The standard assumption 
of the equivalence of relaxation when the magnetic field ro- 
tates about a grain or a grain rotates in a static magnetic 
field is incorrect; the difference is directly related to the 
spontaneous magnetization due to the Barnett effect. 

Although present for all grains, resonance relaxation is 
most prominent for the smallest ones. When grains rotate 
very rapidly, as is the case for very small grains, the reso- 
nance relaxation effect ensures that \" does not plunge as 
the rotation frequency increases. As a result, we conclude 
that small grains (e.g. a < 10 cm) should be param- 
agnetically aligned. The degree of their aligment depends 
on the particular phase of the interstellar medium and on 
the efficiency of spin-lattice relaxation. The latter factor 
is unfortunately uncertain for very small grains for which 
the existing laboratory data is not applicable. 

If the ultrasmall grains are partially aligned, the implica- 
tions arc as follows: (1) The microwave radiation described 
in DL98ab will be polarized - by a few % - and could have 
dramatic consequences for experiments - such as MAP or 
PLANCK - designed to measure polarization of the cos- 
mic microwave background. (2) If the grain body axes are 
aligned with J, then absorption by these small grains will 
contribute to starlight polarization in the ultraviolet, and 
(3) the infrared emission following absorption of starlight 
photons by these small grains will also be polarized. How- 
ever, the contribution to starlight polarization is expected 
to be small due to only partial alignment of the grain body 
axes with J. The infrared emission will be even less polar- 
ized, due to disorientation of the grain axes (Lazarian & 
Robcrgc 1997) during the thermal spike following a pho- 
ton absorption, i.e., while the infrared emission is taking 
place. 
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Fig. 1. — Measure of grain alignment for both resonance relaxation and Davis-Grccnstcin relaxation for grains in the cold interstellar 
medium as a function of frequency (left) and size (right). For resonance relaxation the saturation effects (see eq. 8) are neglected, which 
means that the upper curves correspond to the maximal values allowed by the paramagnetic mechanism. The amplitude of discontinuity on 
the left is smaller than that on the right as both spherical and planar grains contribute to the microwave cmissivity at 20-30 GHz. 



